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Abstract 

Ql I For Q G (1, 2) we consider the equation dtu = A" "m — rb ■ Vu, where 

CLj ' 6 is a divergence free singular vector field not necessarily belonging to the 

Kato class. We show that for sufficiently small r > the fundamental 
solution is globally in time comparable with the density of the isotropic 
Cy ■ stable process. 
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^r^ ! 1 Introduction 



en , 

I • I Let (i > 1 be a natural number and a S (1, 2). We denote by p(t, x) the density 

;<— ^ ■ of the isotropic a-stable Levy process, i.e. 

p(i,x) = (27r)-'^ / e-"-«e-*'«'"de, t>0,xeR'^. (1) 

K> , For (f e C^(R ) we define the operator 

$^ ' 

A / </)(a;) = ^d,a hm / — — ^ dy , 

where Ad,a > is a constant depending only on a and d. A"'^ is the in- 
finitesimal generator of the isotropic a-stable process with the time and space 
homogeneous transition density p{t, x, y) = p{t, y — x), i.e. 

A"/20(x) = lim i / pit,x,y){^iy)-^ix))dy. 
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Let b{x) — {bi{x), . . . , bd{x)) be a vector field satisfying tlie following condi- 
tions 

SUV [ p{t,x,y)\b{y)\dy<Cbe^"-\ iG(0,c»), (2) 

divb = 0, (in the sense of distribution theory), (3) 

for some constant Cb > 0. We will study the equation 

dtu- A°'/^u + rb-Vu^O, xeR'^,t>0. (4) 

The condition 1^ allows for functions b not belonging to the usual Kato class 
^rf~^' ^®^ ([7] below). For example our results apply to d = 2 and 

%) = (y2i2/r",-2/ii2/r"). 

We note that div& = in the sense of distributions, \b{y)\ = ly]^^" and b ^ 
^d~^- "^^^ main result of the paper is the following. 

Theorem 1. There is a constant R — R{a, d, Cb) such that for all < r < R 
there exists a function p{t, x, y) such that for (j) G C^(R, M.'^), s € R and x £ M.'^ , 

p{u—s, X, z){du(j){u, z)+A" "(/i(u, z)+r6(z)-Vz(/)(u, z))dzdu — —(j){s^ x) , 

(5) 
and there is a constant C > depending only on d, a, 6, R such that 

C-^p{t,x,y)<p{t,x,y)<Cp{t,x,y), t>{),x,y&W'. (6) 

According to dS]), p is the integral kernel of the left inverse of —{dt + A"' + 
rb ■ Vz). Put differently, since div6 = 0, a function /: {u,z) i— )■ p{u — s,x,z) 

solves (dt — A"' + rb ■ V z)f = 5{s.x) in the sense of distributions. Thus, p is the 
fundamental solution of (|3]). As a corollary we obtain that p is the integral kernel 
of the Markov semigroup with the (weak) generator A"/^ + & • V ( Corollary [T3]). 

Equations similar to ^ were widely studied for the Laplacian and more 
general elliptic operators (see e.g. [2], [19], [20], [H]). The authors considered 
also drifts b depending on time and satisfying various conditions similar to 
([T]) and called the generalized Kato conditions . In general, the comparability 
of the fundamental solution with the Gaussian density holds only locally in 
time (i.e. with constants depending on time). To obtain global estimates, the 
additional assumption on divergence of b is necessary. For example in [17] Osada 
proved that the fundamental solution of dtu — Au + b ■ Vu has upper and lower 
Gaussian bounds, where A is a uniformly elliptic operator in divergence form, 6 
is the derivative of a bounded function and div6 = 0. These results were later 
obtained in |15| for much more singular drifts (exceeding generalized Kato class) 
with some smallness assumption on divergence. 

Additive perturbations of the fractional Laplacian were intensely studied 
in recent years (e.g. [5], [ID], [i], [13], [1], [7], [5], [T], [5]). In particular. 



the equation Q was considered in [5] for b G K.^ but with no condition on 
divergence. Recall that b S IC^~^ if 

lim sup / / s^^^°'p{s,x,y)\b{y)\dyds ^0 (7) 



(see also [12] and [13] for further developments). The authors obtained local in 
time comparability of p and p for each r G M. The function p was constructed 
as the perturbation series p — X^^o^"' where 

Po{t,x,y) ^p{t,x,y), (8) 



Pn{t,x,y)^ I I pn-i{t - s,x,z)b{z) ■V;,p{s,z,y)dzds. (9) 

Jo JRi' 

The Kato conditions on b assures local in time smallness of pi with respect to 
p, and in consequence the perturbation series converges. Similar methods was 
used to study Schrodinger perturbations of transition densities ([3], [TT]) and 
the Green function of A"/^ + b{x) -V^ ([4]). 

In the present paper we will also use the technique of perturbation series, but 
in our case the conditions on h only assures the finiteness of pi (see Lemma [5]) . 
In addition, the integral in Q is not absolutely convergent as the integral over 
time and space, which makes the proofs much more complicated and delicate. 
In particular, it is not obvious that the functions p„ are well defined. In order 
to prove it we introduce functions Pnit,x,y,s,z) (see ([20)1 for definition) and 
integrate it separately over (M'')" and n-dimensional simplex S'„(0,i). Namely, 
we consider 



\p\nit,x,y) = 

'S„(0,t) 



Pnit,x,y,s,z)dz 



ds. 



the majorants of the functions p„ (see (l22t and (|35|) ). To estimate |p|„ for 
n > 2 we split the integral over the simplex S'„(0,i) into suitable n + 1 parts. 
As a consequence Motzkin numbers appear in the estimates of p„ (see [11] for 
another connection of the perturbation series with combinatorics). In order to 
assure the convergence of the perturbation series, smallness of pi is needed and 
this is why we multiply b by small constant in Theorem [TJ We like to note that 
Theorem [1] should hold for i? = oo, but such an extension calls for different 
methods. 

One of the tool used in this paper is so called 3P theorem (see [5], [12], [T3]). 
It allows to suitably split a ratio of three functions p, and in consequence to 
estimate p„. Since for a = 2 (Gaussian case) 3P theorem does not hold, our 
method cannot be applied to perturbations of the classical Laplacian. Similarly 
as in previous papers we also exclude the case of a < 1. Although Lemma [3] 
holds in this case. Lemma H] does not seem to extend to a < 1 and consequently 
our approach does not work for this case. 

The paper is organized as follows. In Section 2 we collect basic properties of 
the transition density p{t,x,y). In Section 3 we define and estimate functions 
Pn- In Section 4 we prove Theorem [TJ 



All the functions considered in the sequel are Borel measurable. When we 
write f{x) « g{x), we mean that there is a number < C < oo independent 
of X, i.e. a constant, such that for every x we have C^^f{x) < g{x) < Cf{x). 
As usual we write a Ab = min(a, 6) and ay b — max(a, 6). The notation 
C = C{a, b, . . . ,c) means that C is a constant which depends only on a, &,..., c. 

2 Preliminaries 

Throughout the paper d > 1, and unless stated otherwise, a € (1, 2). In Lemmas 
[11 [21 13] we recall well-known results about the density p{t, x, y) of the isotropic 
d-dimensional a-stable process (see [5] for details). 

Lemma 1. There exist a constant Ci such that 

t 



CI' 



^-d/a ^ 



.\d+a 



<p{t,x)<Ci 



-d/a 



.\d+a 



te (0,cx)), xeR'^. 
(10) 



Lemma 2 (3P). There exist a constant C2 such that 

p{t,x,z)p{s,z,y) < C2p{t+s,x,y)[p{t,x,z)+p{s,z,y)], s,t e (0, cx)), x,y,z e M'' 

Let p^™"' be the a-stable density in dimension m. 

Lemma 3. For all t > and x £ M."^, 

V^pid) {t, x) = -2^xp('^+2) (<, x) , (11) 

where x G M'^^^ is such that \x\ = \x\. 

Applying dTU]) to ^Q we get 

l^xPit, x)\ < CT,t-^l°'p{t, x) , i e (0, 00), X € R'' (12) 

Our aim is to prove that functions p„ defined in ([8]) and ((9]) satisfy pn{t,x,y) < 
Cnp{t,x,y), where C„ are the constants such that J2^=oCn < co. Since the 
condition © does not guarantee convergence of the integral 

/ p{t~s,x,z)\b{z)\\'Vzp{s,z,y)\dzds, (13) 

we cannot follow the proofs from [S] and [13] . However the inner integral of ([1 
is convergent. Indeed by p^ . Lemma [Hand ([5]) we have 



p{t~ s,x,z)\b{z)\\\'zp{s,z,y)\dz 

< CiS~^'" / p{t ~ s,x,z)p{s,z,y)\b{z)\dz 

<C2S~^^'^p{t,x,y) / {p{t- s,x,z)+p{s,z,y))\b{z)\dz 

< C2CfcS-i/"[(t - s)i/"-i + si/"-i]p(t, X, y) (14) 



Therefore instead of (IT^ we will consider 



p{t — s,x,z)b{z) -V zP{s.,z,y)dz ds . 

In order to estimate (1151) we will use the following lemma. 
Lemma 4. For all s,t > and x,y eR"^ we have 



(15) 



p{t,x,z)b{z) ■\'zPis,z,y)dz ^ - I \'zp{t,x,z) ■ b{z)p{s, z,y)dz . (16) 
Proof. Let g E C;?°(M'*) be such that 



9{z) = 



1 for \z\ < 1, 
forU|>2. 



Let /„(z) = g{z/n)p{t, x, z)p{s, z, y) e Cf^(R'*). Then V Jn{z) ^ ^ ^{p{t, x, z)p{s, z, y)) 
as n -> oo. Furthermore 

|V^/„(z)| < c(|V3(p(i,a;,z)p(s,z,y))| + p{t,x,z)p{s,z,y)) 

for some constant c > 0. By p^ . Lemma [5] and (I2|) 



[|V^(p(i,x,z)p(s, z,j/))| +p(t,x,z)p(s,z,y)]|6(z)|(i2; < oo . 
Therefore by ([3]) and Lebesgue theorem 

= lim / Vzjn{z) ■b{z)dz ^ / 'Vz{p{t,x,z)p{s,z,y)) ■ b{z) dz , 



which ends the proof. 

Similarly condition ([3]) yields that for all s,t > and ^,y G 

[&(0 ■ V^p(t, e, ^)]K^) • V.p(s, z, y) dz 

V, [6(0 • V5p(i, e, ^)] • b{z)p{s, z, y) dz . 



U 



(17) 



In the following lemma we will use (1161) to show that the function pi introduced 
in (J9l) is well defined. In a similar way we will apply (J17p to estimate other 
functions p„. 

Lemma 5. There exists a constant C such that for all t > 0, x,y £ M."^ , 



p(t - s, X, z)b{z) ■ VzPis, z, y) dz 



ds < Cp{t, X, y) . 



(18) 



Proof. By Lemma S] and ([HI) we obtain 



p{t - s, X, z)b{z) ■ Vzp{s, z, y) dz 



ds 



tji 



t/2 
tji 







t/2 



p(t — s, X, z)6(z) • Vzp(s, z, y) dz 

d 

p{t - s, X, z)5(z) • V2p(s, z, y) dz 
V^p(t - s, x, z) • b{z)p{s, z, y) dz 

d 

p{t - s, X, z)5(z) • VzP{s, z, y) dz 



ds 
ds 

ds 
ds < cp[t,x,y) 



D 



We will also need the following two auxiliary lemmas 



Lemma 6. There exists a constant C4 such that for all t > and z,w (£ M.'^ we 
have 

\b{z) ■ WMw) ■ V^p(t, z, w))\ < Ci\b(z)\\b{w)\p{t, z, w)t-^'" . 

Proof. From (fTTj) we get 

V,„p(''' (t, z, u;) = 27r(z - w)p^'^+'^'^ {t,z,w), 
V,p^'^+^\t, i, w) = -2tt{z - w)/^+^\t, z, w), 



where z ~ (z, 0, 0) G R'^+^ and z = (z, 0, 0) G R'^+'' (we use similar notation for 
It; and w). Therefore, 

b{z) ■ Wz{b{w) ■ W^p{t, z, w)) = 27r6(z) • \/,[b{w) ■ (z - w)p('^+^'> {t, z, w)] 
= 27r6(z) • [6(w)p(''+2)(i,z» - 27r(6(w) • (z - u;))(z - w)p'^'^+^\t, z,w)] 
= 27r6(z) • 6(w)p(''+2)(t, i, i&) - 47r2(5H • (z - w))6(z) • (z - w)p^'^+'^'>{t, z, w) 

Applying pO|) we obtain the assertion of the lemma. D 

Lemma 7. There exists a constant C5 such that for all t > 

t/2 pt 

(r - u)-i{{t - r) = -1 + (r - u) = "^)((t - ii)°"^ + u°"^) drdu < C5 . 



Jt/2 



Proo/. First we note that aP + bP < 2^^P{a + b)P for a, 6 > and < p < 1. 
Consequently a~P + b^P < 2^^P{a + b)P{ab)~P. Hence it suffices to show 



t/2 ft 



Jt/2 



{r - u) " ^(t-r)" ^u" ^t^ "drduKc. 



Splitting the second integral into intervals (t/2, 3i/4) and (3i/4, t) we get 

i-t/2 Pt ^ ^ ^ 1-1/2 

Jo Jt/2 Jo 

D 

We note that Lemma [5] extends to a < 1, however the following lemma does 
not and this is why we generally assume a £ (1,2) in the paper. Lemma |S] will 
allow us to estimate the functions Pn for n > 2. 

Lemma 8. There exists a constant C such that for all t > and x,y £ R*^, 

ft/2 pt 



dw dS, dr du 

lo Jt/2.' 

p{u,x,(,)\W^{b{0 • V5p(r-u,^,u;)) ■ b{w)\p{t - r,w,y) < Cp{t,x,y) . (19) 
Proof. By Lemma [2] and ^ 

p{s,x,z)\b{z)\p{t,z,v)dz<cp{s + t,x,y){s^l^'^+t^'''~^). 

Hence by Lemma [6] we have 

p{u, X, ^) I V^^, {h{i) ■ V^p{r - u, ^, w)) ■ h{w) \p{t - r, w, y) dw d^ 

< I p{u,x,£.)\b{0\{r^u)~^/°'p{r-u,^,w)\b{w)\p{t-r,w,y)dwd^ 

JW Js^d 

< I \b{£.)\{r~u)~'^/°'p{u,x,(,)p{t~u,i,y){{t-rY'°'-^ + {r-uf/°'-^)dw 

< p{t,x,y){r - M)-2/"((t _ r)i/"-i + (r - w)i/"-i)((t ~ u)^^"-^ + u^/""^) . 
Now ([T9l) follows from Lemma [71 D 

3 Perturbation series 

In this section we introduce functions |p|„ which will be majorants of the func- 
tions Pn (see ^). For any a < b and n > 1 we denote 

S„{a, b) = {(si, S2, . . . , s„) e M" : a < si < S2 < . . . < s„ < &} . 

For any i > and x, y € M'' let 

P„(i, X, y, s,z) = p{si,x, Zi)b{zi)-V ziP{s2-Si, Zi, Z2) . . . b{Zn)-W^^p{t-Sn, Zn,y), 

(20) 
where s = (si,...,s„) € S'„(0,i) and z = (zi,...,z„) e (R'')" . We recall 
that the integrals in ^ may not be absolutely convergent as the integrals over 
[0,t] X M"*, and the approach from the paper [5] cannot be used. Therefore we 
separate integrals over time and space in the following definition. 



Definition 9. For any f > and x,y eM.'^ we define 
\p\o{t,x,y) =p{t,x,y) 



\p\nit,x,y) 
We note that by Lemma H @ and p^ 



Pnit,x,y,s,z)dz 



ds. 



/ \Pnit,x,y, s, z)\dz < oo, 



(21) 
(22) 

(23) 



hence the functions |p|„ are well-defined (at most they are equal to infinity). 
However the integral /^ /g ^x /md\n !-?«(*, x, y,s,z_)\ dzds may not be convergent 
because singularities of the gradient of the functions p in pO|) may be not inte- 
grable in the whole simplex S'„(0, t). Therefore in order to estimate ([22]) we use 
the following representation 



/'n-l 



SniO,t) = Sn{0,t/2)U U Sn^ki0,t/2) X Skit/2, t)\ USnit/2,t), (24) 

Lemma |4] and (1171) to move these singularities off the region of integration. 
Lemma 10. For any l<k<n — l,t>0, and x,y E M."^ we have 

Pnit,x,y,s,z_) dz_ ds 

/S„_fc(0,t/2)xSfc(t/2,t) 

< f [ [ f dwd^ 

Jo Jt/2 Jk<* Jr'^ 

\p\n-k-i{u, X, 0\^w (HO ■ ^iP{r - u, C, w)) ■ b{w) I \p\k-i{t - r, w, y) dr du 



(25) 



Proof. By (|23p and Fubini's theorem we may change the order of integration in 
integrals over (R'')". We note that 



\p\rnit-r,x,y) 



Prnit-r,x,y,s,z)dz 



ds . 



(26) 



Changing the order of variables in the following way, 



f{s)dSmdSm~l .. -dsi = / f{s)dSm^l . . . dsidSr. 

Srr,{a,b) J{a,b)xS^-i(a,s^) 



using (IT71) and Fubini's theorem we get 



/ 



S„_fc(0,t/2)xSfc(t/2.t) 
t/2 






ds 



duduk drdrn_f; 

JS„_fc_i(0.r„_fc) Jt/2 JSfc_i(«fc,t) 

6(wi) • Vu,iPfc_i(t - ui,wi,y,u,w)dwdwidzdzn-k 

I II dudukdrdrn^k 

JS„_fc_.i(0,r„_fc) Ji/2"'S'fc_i(«fc,t) 

Pn-k-lirn-k,X,Zn-k,L,z)Pk-l{t - Ui,Wi,y,U,w) 



b{wi) ■ Vnii [b{zn-k) ■ Vz„_^p(ui - r„_fc, z„_fc, wi)] dwdwi dzdzn-k 



d^n—k—l 

n 



where r = (''i,- • -j^n-fc-i), u = (mi,. . . jW^^i), z = (zi,. . .,z„_fe_i), w 
{w2, ■ ■ ■ , Wfc). Now sphtting mtegral over (R'')" into integrals over (R'^) 
(E'*)'=-i, and (M'')2 and applying ^, we get (gl]) 

Lemma 11. For n > 2, t > and x,y dR"^ we have 

\p\n{t,x,y)< / \p\n-i{u,x,z)\b{z)-Vzp{t-u,z,y)\dzdu 

Jo JR'' 

+ / I'^zPiu^x^z) ■b{z)\\p\,i^i{t-u,z,y)dzdu 

J t/2 JM'' 

"-2 „t/2 ft r f 
+ ^ / III dw dzdr du\p\k{u,x,z)x 
fe=0 " "^'/^ •^*'' •^*'' 

X |6(z) • V2(6(it;) • VwP{r- u, z,w))\\p\n-2-k{t - r,w,y) . 

Proof. By ^^ we get 



b|n(i,a;,2/) 



S„(0,t/2) 
S„(t/2,t) 



/ Pnit,x,y,s,^)d^ 
/ Pn{t,x,y,s,^)d^ 

J(M'')" 



ds 
ds 



+ 



UkZl S„-k{0,t/2)xSkit/2,t) 



/ Pn{t,x,y,s,Qdi 

J(R'')" 



(27) 

(28) 

ds (29) 



The integral (l?71) is estimated as follows, 



/ / Pnit,x,y,s,C)d§^ 

JS„(0.t/2) J(R'')^ 



'S„(0,t/2) 
f.t/2 



"'S„_i(0,s„) 
t/2 



Jk<* 



Pn-l{Sn, X, £,n, S* , CMCn) ' ^ i„P{t - Sn,£,,i,y)d^ 

df „ ds* ds„ X 

J(R<i)"-i 

|p|„_i(s„,a;,^„)|6(^„) • V^^p{t- Sn,^n,y)\d^ndsn 



ds 



< 

/O "'S„_i(0,s„) 

X l^(C«)Ve„p(i-s„,C„,2/)| 

y.t/2 

< 



where s* = (si, . . . Sn-i) and ^* = (^i, . . . Cn-i)- Applying Lemma |4] and using 
similar method we estimate ([28| . Next, we split ([29|) into n — 1 integrals over 
the sets 5„_fc(0,i/2) x Sk{t/2,t) and apply Lemma [TUl to each integral. D 

By the lemmas from the previous section and induction we will obtain that 
all functions \p\n are finite and in consequence the functions p„ are well defined. 
Detailed estimates will we given in the next section. 



4 Proof of Theorem 1 

Before we pass to the proofs of the main theorem we briefly introduce the 
Motzkin numbers. In combinatorics Motzkin number M„ represents the number 
of different ways of drawing non-intersecting chords on a circle between n points 
([16]). Their generating function is (see [18]) 



M{x) = Y^ MnX" = 



„ 1 — X — \/l — 2a; — 3a 



n=0 

and the following recurrence relation holds, 

n-2 



2a;2 



(30) 



M„ = M„_i + ^MfeM„_2-fe, Mo = Afi = l 



(31) 



fc=0 



Wc may now prove the main estimates of this paper 
Lemma 12. There is a constant C such that for all t > Q, x,y E M. and n > 1, 

\pUt,x,y)<M^Cy{t,x,y). (32) 



10 



Proof. Let ci be the constant such that (see the proof of Lemma [5|) 

f.t/2 f 

p{s, X, z)\b{z) ■ V zp{t — s, z, y)\ dz ds 
+ / / \'^ zp{s,x,z) ■h{z)\p{t~ s,z,y)dzds <cip{t,x,y) . (33) 

Jt/2 Jm'' 

Let C = ci V y/c2, where C2 is the constant from Lemma [8l 

We use induction. For n = 1 we apply Lemma O Suppose ([32l) holds for 
n = 1, . . . fc — 1. By Lemma [TT] we get 

\p\k{t,x,y)<C^^^Mk-i i I p{u,x,z)\b{z)-V,p{t~u,z,y)\dzdu 

Jo JR<* 

+ C'''^Mk-i / \S/zp{u,x,z)-b{z)\p{t~u,z,y)dzdu 

Jt/2 Jw 
k-2 ^t/2 ft 

• du 
p{u,x, z)\b{z) ■ S/zibiw) ■ S/wP{r ~ UyZ^w))] x p{t — r, w, y) 



y^C^MjC^^^~^Mk-2-j / / dwdzdri 

-_n Jo Jt/2 JR'' JRt' 



Now by (|33|) , Lemma |8] and ([31]) we obtain 

\p\k{t, x,y)<l C'-'ciMk-i + J2 C=M^C^-^-^Mk-2-jC2 j p(i, a^, v) 

< C'' I Affc_i + J2 MjMk-2-j I p(i, X, y) = C''Mkp{t, x, y) . 

U 

For all n £ N let p„ be functions satisfying (j8]) and ([9]) . 

Corollary 13. Functions p„ are we^/ defined and there is a constant C such 
that for all t > Q, x,y G M.'^ and n > 1 

\pn{t,x,y)\<M„Cy{t,x,y). (34) 

Proof. We simultaneously prove the estimates of Pn and that they are well 
defined. It suffices to show that for n > 1 

Pn{t,x,y)= / Pn{t,x,y,s,z)dzds, i > 0, a;, y G K"* . (35) 

JS„(0,t) J(R<i)" 

We use induction. For n — 1 psp matches the definition of pi. Suppose psp 



11 



holds for n eN. By Lemmas [T^ [^ and [51 be have 



/ / / Pnit- u,x,£_,s,z)dz 



mmiPi^,^,y)\dsd^ 



JR'' 

< c{{t ~ u)i/"-i + wi/"-i)u-i/>(i, X, y). 
Therefore by Fubini's theorem and (155)) 
Pn+i{t,x,y) 

Pn{t - u, X, Cjb{Cj ■ V5p(u, f , y)d^du 



(36) 



JR<* 



III I Pn{t- u,x,^,s,z)dzdsb{^) ■^^p{u,£_,y)d£_du 

Jo Jr-^ Js„{0,t-u) JiR-^)" 

/ / / / Pnit~u,x,^,s,z)b{^)-V^p{u,^,y)dzd£,dsdu 

Jo Js„(0,t-u) Jr'^ J(R'')" 



,(0,t-M) 

<\p\n+i{t,x,y), 
which ends the proof 



D 



Proof of Theorem d Let C be the constant from Corollary [131 Let r be such 
that rC < (V^-l)/4. Then 



\pnit,x,y)\ < A'UrCrpit,x,y) ^ 
Denote r] — rC. We define p{t, x, y) as 

p{t,x,y) = ^pn{t,x,y) 

n=0 

By ([57)1 and (1501) . the series converges, and 



(37) 



(38) 



l-f]- Jl-2n- irf , 
P{t, X, y) < — p(t, X, y) . 



Furthermore, 



Pit, X, y) > p{t, x,y)-^ |p|„(t, X, y) > ^^-:^ p{t, x, y) 



n=l 



2r?^ 



We next prove that for e C;?°(R, R''), s e R and x e R'*, 

p{u — s,x, z){du4>{u, z) + A""0(u, z) + 6(z)- V20(u, z))dzdu — —(j){s,x). 
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By dnHI we get 

OO 

pit, X, y) = p{t, x,y) + ^ pn{t, X, y) 

OO „f „ 

^p{t,x,y) + '^ / pn-i{t- s,x,z)b{z) ■\/zp{s,z,y)dzds 

= p{t,x,y)+ / p{t~ s,x,z)b{z)-Vzp{s,z,y)dzds (39) 



Here Fubini theorem is justified by similar arguments as in the proof of (|35|) . 
The rest of the proof is the same as in [T31 Theorem 1] D 

Corollary 14. The function p satisfies the Chapman-Kolmogorov equation 

p{s,x,z)p{t,z,y) dz = pit + s,x,y), s,t>0, x,y£M. , (40) 

and a family of operators 

Ptfix)^ f p{t,x,y)fiy)dy 



forms a Markov semigroup with a (weak) generator A"'^ + b{x) ■ Wx- 

Proof. For the proof of (gO]) see [H Lemmas 15, 16]. By ((391), dH) and Fubini 
theorem 

p{t,x,y)dy^ ip{t,x,y)+ / p{t - s,x, z)b{z) ■ \7 ^p{s, z,y)dzds) dy 

I JR<i V Ja Jr'' J 

= 1+/ / p{t- s,x,z)b(z) -^^i I p{s,z,y)dy]dzds^l. 
Now, let /,g e C^{K.<^). We will show that 

lim / ^^/(^) - /(^) ^(^) ^^ ^ /■ (A"/2/(x) + 6(x) • ^f{x))g{x)dx (41) 



5(2;) da; 

g{x) dy dx 



By (1391) 

A/(x) - fix) 



t 
Pit,x,y)ifiy)~fix)) 



^ r r rt 



t JRd JRti Jo 



t 

Pit - s, x, z)biz) ■ ^zPis, z, y)fiy)gix) dz ds dy dx 
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The first summand converges to /^^ A°'/'^f{x)g{x)dx. By careful use of the 
Fubini theorem 



hit) ^ - / / / p{t~ s,x,z)p{s,z,y)b{z)-Vyf{y)g{x)dzdsdydx. 

If we denote by p'^^' the function p perturbed by b then p(i, x, y) = p^^' {t, x, y) = 
p^~''\t,y,x). Hence J^ap{t,x,y)dx = 1 for i > and y G R''. Therefore by (O 



hit)- / biz)-Vfiz)giz)dz 



pit- s,x,z)pis,z,y) 



jR-i JR'i jR-i Jo '• 

X \^yfiy)9ix) - ^ zfiz)giz)\ dz ds dy dx, 



and the last expression converges to as t — >■ (for details see the proof of [3 
Theoreml]). D 
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